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Abstract 



We consider a one-parameter family of third order nonlinear recurrence relations. 
Each member of this family satisfies the singularity confinement test, has a conserved 
■ quantity, and moreover has the Laurent property: all of the iterates are Laurent polyno- 

c~| ■ mials in the initial data. However, we show that these recurrences are not Diophantine 

integrable according to the definition proposed by Halburd (2005 J. Phys. A: Math. 
Gen. 38 LI). Explicit bounds on the asymptotic growth of the heights of iterates are 
obtained for a special choice of initial data. As a by-product of our analysis, infinitely 
many solutions are found for a certain family of Diophantine equations, studied by 
' Mordell, that includes Markoff's equation. 

^1- ' 

For some time there has been considerable interest in maps or discrete equations which 
are integrable. Various different criteria have been proposed as tests for integrability in 
the discrete setting. One of the earliest proposals was the singularity confinement test of 
Grammaticos, Ramani and Papageorgiou ^2], which has proved to be an extremely useful 
tool for isolating discrete Painleve equations (see e.g. [2B])- However, after Hietarinta 
and Viallet's discovery of some non-integrable equations with the singularity confinement 
property, they were led to introduce the zero algebraic entropy condition for integrability 
of rational maps [13] , namely that the degree d n of the nth iterate of a map (as a rational 
function of the initial data) should satisfy lim n _ >00 (log<i n )/ri = 0. The phenomenon of weak 
degree growth for integrable maps had been studied earlier by Veselov [S3], and the 
notion of algebraic entropy proposed in ^3] is connected with other measures of growth such 
as Arnold complexity [1. 

Since the work of Okamoto [22] it has been known that the Backhand transformations 
of differential Painleve equations can be classified in terms of affine Weyl groups. Noumi 
and Yamada have shown that the symmetries of these ODEs are simultaneously compatible 
with associated discrete Painleve equations [2Ij, and this viewpoint has been explained geo- 
metrically by means of Cremona group actions on rational surfaces j3U] • Discrete integrable 
systems are also related to the notions of discrete geometry ||j and discrete analytic functions 


In the past few years, several new techniques have been developed for isolating integrable 
discrete equations: analytical criteria on the one hand, and arithmetical ones on the other. 
The Painleve property has been used very effectively to isolate integrable differential equa- 
tions 7\, and this led Ablowitz, Halburd and Herbst to extend this property to difference 
equations using Nevanlinna theory [2], by considering the asymptotic growth of meromor- 
phic solutions at infinity. Roberts and Vivaldi [2H] have studied the distribution of orbit 
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lengths in rational maps reduced to finite fields F p for different primes p, in order to identify 
integrable cases of such maps. Most recently [T3] . Halburd used Vojta's dictionary between 
Nevanlinna theory and Diophantine approximation to translate the concepts of [2] into a 
Diophantine integrability criterion for discrete equations. 

There is by now a considerable literature on discrete bilinear equations, including the 
bilinear forms of discrete Painleve equations [28J, and bilinear partial difference equations 
such as the discrete Hirota equation, which can be written (with two arbitrary parameters 
ct,/3) in the form 

+ (3 

The partial difference equation (|T} has continuum limits to all of the bilinear equations in the 
KP hierarchy of PDEs J3HJ, and also appears as an equation for transfer matrices in quantum 
integrable models ^H] ■ However, despite the huge amount of theory that has been developed 
for bilinear discrete equations, there is another property of such equations which appears to 
have been overlooked by the integrable systems community, namely the following: for initial 
data specified on a suitable subset of Z 3 , all of the iterates of (JIJ are Laurent polynomials 
in the initial data with coefficients in Z[a, /?]. This Laurent property was originally known 
only to a few algebraic combinatorialists |25., and for the equation (JIJ it was first proved 
by Fomin and Zelevinsky within the framework of the theory of cluster algebras [ID] , while 
more detailed properties of the associated Laurent polynomials have been shown by Speyer 

EH 

One of the simplest manifestations of this Laurent phenomenon was found by Michael 
Somos, who considered bilinear recurrences of the form 

[fc/2] 

T n+k T n = ^ T n+k-j T n+j, k > 4 (2) 

taking the initial values To = T\ — . . . — t^-i = 1. Clearly for the recurrences (J2J) each new 
iterate r n+ k is a rational function of the initial data, so with this particular choice one would 
expect r n to be rational numbers, but it was observed numerically that for the Somos-4 
recurrence 

T n +4 T n — Oi Tn+3 r n+l + /3(r n+2 ) 2 (3) 
with parameters a = (3 = 1 and starting with four ones, an integer sequence 

1, 1, 1, 1, 2, 3, 7, 23, 59, 314, 1529, 8209, 83313, ... (4) 

results. Similar empirical observations showed that Somos-5, -6 and -7 also yield integer 
sequences, while the Somos-fc recurrences for k > 8 do not. The first inductive proof of 
the integrality of the sequence (JD) appeared in [XT, but it was realized that the deeper 
reason behind this lay in the fact that the recurrence (0) has the Laurent property, that is 
r n £ Z[q;, /3, Tq 1 , r-f 1 , T2 1 , Tg 1 ] for all n E 7h. In other words, the iterates of © are Laurent 
polynomials in the four initial data whose coefficients are in Z[a, so if a — (3 — To = T\ — 
t~2 = T3 = 1 then the sequence must consist entirely of integers. 

Fomin and Zelevinsky found a suitable setting within which to prove the Laurent property 
for a variety of discrete equations [ID] , including certain recurrences of the form 

T n+k T n = F(r n+ i, . . . , T n+ fc_i 

) (5) 

for particular polynomials F (mostly, but not all, quadratic forms in their arguments), as 
well as integrable two- and three-dimensional bilinear recurrences like (0). More recently, the 



connection between the iterates of Somos-4 recurrences and sequences of points on elliptic 
curves was explained in the PhD thesis of Swart [H2J, while independently the author found 
the explicit solution of the initial value problem for both Somos-4 ^H] and Somos-5 ^Bj 
in terms of elliptic sigma functions; see also [21] for a connection with continued fractions. 
An essential observation in ^S] and JHj was that each of these bilinear recurrences (fourth 
and fifth order respectively) could be understood in terms of a suitable integrable mapping 
of the plane, with each one corresponding to a particular degenerate case of the family of 
maps studied by Quispel, Roberts and Thompson j2Zj. Somos sequences are of considerable 
interest to number theorists due to the way that new prime factors appear therein |Hl E] • 

A natural question raised by all of the above is the following: if discrete equations with 
the Laurent property like (JTJ) and @ are integrable, then are all discrete equations with the 
Laurent property integrable? The converse is clearly not possible: most integrable rational 
maps do not have iterates that are Laurent polynomials in the initial data. Thus it is sensible 
to restrict ourselves to a suitable sub-class of maps, namely those defined by recurrences of 
the form The purpose of this note is to show that the answer to this question is a 
negative one, by taking an example not considered in JU] > namely the one-parameter family 
of third order recurrences given by 

T n +3T n = Tl +2 + Tl +1 + J. (6) 

In the case J = 0, it was noticed by Dana Scott [11] that with r = T\ = r 2 = 1 an integer 
sequence results: 

1,1,1, 2, 5, 29, 433, 37666, 48928105, 5528778008357, ... (7) 

In fact, © is related to some old problems in number theory: the numbers (J7J) are a sub- 
sequence of the Markoff numbers that arise in the theory of indefinite quadratic forms [T9] . 
while the case of arbitrary J is related to a Diophantine equation considered by Mordell j2D3 • 
Below we give two different proofs that the recurrence © has the Laurent property; the first 
proof also shows that it passes the singularity confinement test, while the second one gives 
a stronger result based on the existence of a conserved quantity. However, when Halburd's 
Diophantine integrability test is applied to a particular class of sequences generated by ©, 
it is seen that the growth of the heights of iterates is double exponential (which can already 
be guessed from the rapid growth of the terms in (J7J)). For the the particular sequences being 
considered, quite sharp bounds on this double exponential growth are obtained. Hence we 
have an interesting example of a non-integrable recurrence with the Laurent property. 
The recurrence © defines a rational map of three-dimensional affine space, i.e. 

(x n ,y n ,z n ) i (y n ,z n , (y 2 n + z 2 n + J)/x n ) = (x n+1 , y n+1 , z n+1 ), (8) 

upon setting (r n , r n+ i, r n+2 ) = (x n , y n , z n ). Before considering the properties of this map 
in more detail, we present the first proof that it has the Laurent property, since this will 
simultaneously show that it satisfies singularity confinement as well. Letting r = a, T\ — b, 
T2 = c and putting the recurrence into MAPLE, we see that the next two iterates are 

b 2 + c 2 + J b 4 + a 2 c 2 + 2b 2 c 2 + c 4 + (a 2 + 2b 2 + 2c 2 ) J + J 2 

n = , n = n-r , 

a a A b 

and T5 has a 4 b 2 c in the denominator and its numerator is a polynomial in Z[a, b, c, J] with 
29 terms. The first miracle occurs when n — 3, because upon dividing the right hand side of 
(0) by r 3 the numerator b 2 + c 2 + J cancels, so that r 6 G 1Z = Z[a ±:L , 6 ±:L , c ±:L , J] with a 7 6 4 c 2 
in the denominator and a polynomial with 104 terms in the numerator. 



Now to prove the Laurent property, we can use the fact that modulo monomial factors 
a?b k c the ring 7Z of Laurent polynomials is a unique factorization domain (see JT0] - the case 
J = is covered by Theorem 1.10). Thus we can consider divisibility of the r n modulo such 
monomials. As the inductive hypothesis, assume that Tj G 1Z for < j < n + 5, that any 
three adjacent terms Tj,Tj + i,Tj + 2 are pairwise coprime (modulo monomials) in this range, 
and also that t 2 ±1 + J is coprime to Tj. To prove that the next iterate r n+6 is in 7Z, consider 
the three adjacent Laurent polynomials x = r n +i, y = r n+2 , z = r n+3 . From successive 
applications of (|BJ) we have 

x 2 + y 2 + J = r n+3 T n = mod z, (9) 

as well as r n+4 = (y 2 + J) fx and r n+4 = (y 4 + J{x 2 + 2y 2 ) + J 2 )/(x 2 y) (both modz), which 
together imply that 

(y 2 + J){x 2 + y 2 + J){y 4 + J{x 2 + 2y 2 ) + J 2 ) 

T n +f>T n+3 = — = mod z (10) 

x A y 2 

by ©, and hence r n+6 G 1Z. Next suppose that an irreducible factor p G 1Z is such that 
p\ r n+6 an d p|T n +5; then p|(r n+6 r n+3 — r^ +5 ) = r^ +4 + J, which contradicts the coprimality 
of T n+5 and r^ +4 + J, so r n+6 and r n+5 must be coprime. The same argument shows that 
r n+6 and r n+4 must also be coprime, so that the three adjacent terms r n+4 , r n+5 and r n+6 are 
pairwise coprime. Similarly it follows that t 2 +5 + J is coprime to r n+6 , and also r 2 +6 + J is 
coprime to r n+ 5, which completes the inductive step. This proves that r„ G K for all positive 
values of the index n, and the same conclusion holds for negative indices by the reversibility 
of the recurrence. 

To see why singularity confinement also follows from the preceding argument, observe 
that a singularity can only occur when one of the iterates vanishes. Thus, without loss of 
generality suppose that z = r n+3 vanishes, and we take r n = w and y 2 = —x 2 — J + ew so 
that z = e — > 0. The above expressions can quickly be evaluated in this limit by noting 
that all terms congruent to z are 0(e), and also y = y + 0(e) with y = \/ —x 2 — J. Thus 
Tn+4 = + O(e), r n+5 = (x 2 + J)/yo + 0(e), and clearly from (fTUj) we have r n+6 e = 0(e) 
so r n+ 6 = 0(1) and the singularity is confined. This argument for confinement will work for 
many other recurrences with the Laurent property, such as those in ^Hj • For the recurrence 
(0), a stronger version of the Laurent property actually holds, as a consequence of the 
following 

Proposition 1 The quantity N n = {j 2 +2 + t 2 +1 + t 2 + J) / '(T„,+ 2 T n+1 T n ) is an invariant for 
the recurrence so that N n = N = N with 

abc 

Moreover, for the same initial data r = a, T\ — b, r 2 = c, the iterates of (0J) are the same as 
those of the recurrence 

T~n+3 = Nr n+2 T n+ i - T n , (12) 

and all r n lie in the ring Z[a, b, c, N] C 1Z — Z[a ±:L , b^ 1 , c ±l , J]. 

The proof is quite straightforward. The fact that N n is a conserved quantity for (0) 
follows from the direct calculation that iV„ +1 — N n = 0, so clearly the iterates of the map 
(151) lie on the cubic surface 



x 2 + y 2 + z 2 - Nxyz + J = 0, 



(13) 



which is non-singular for J ^ 0, — 4/iV 2 . This surface admits some obvious symmetries: it is 
unchanged by permutation of (x,y,z), and also for fixed y and z it is a quadratic equation 
in x and hence is preserved by the involution X : (x, y, z) \— > (x\ y, z) where 



x' 



{y 2 + z 2 + J)/x = Nyz - x 



(14) 



is the other root of the quadratic. Taking the cyclic permutation V : (x,y,z) i— > (y,z,x), 
then performing X followed by V yields the transformation (x, y, z) \— > (y, z, a;'). If the first 
equality for in (JTljl is used in this transformation, then the rational map (JBJ) corresponding 
to the recurrence © arises. Similarly, the alternative recurrence (jl2j) comes from the second 
equality in (|14|). and one could start with this and observe that it has the conserved quantity 
J n = NT n+ 2T n+ iT n — t 2 +2 — t 2 +1 — r 2 . It is also obvious that the iterates of (fl~2|) lie in 
Z[a, 6, c, N], and with iV given by the formula (fTTj) this is a subring of 7?.. 

Here we should point out that the aforementioned symmetries were used extensively by 
Mordell [20J when he considered the Diophantine problem of finding integer triples (x, y, z) 
satisfying (fT3*j) for fixed integers N and J. The case iV = —2, J = —n e Z was treated in 
great detail, since this is related to the problem of finding which integers can be written as 
the sum of four cubes. However, he remarked that "/ know of no test except trial for the 
existence of solutions for given n." The particular case N = 3, J = is another distinguished 
instance of ()13|) known as Markoff's equation, which is related to the spectrum of indefinite 
binary quadratic forms [19J. Although a great deal more is known about this case, there are 
still difficult open problems associated with it [HI El 122] • The case J = has also acquired 
a combinatorial interpretation very recently ^HEE]- The above Proposition implies that, 
given any initial solution triple of integers (a, b, c), the recurrence (jOJ (or (|T2*|) . which is 
equivalent) will generate infinitely many such triples provided that it leads to neither a zero 
nor a periodic orbit. In fact, by the above discussion of singularity confinement, a single zero 
(e.g. a triple (x, y, 0)) is not a problem in the sense that one can introduce a small parameter 
e and analytically continue through the apparent singularity. However, a solution like (x, 0, 0) 
can also occur if —J is a perfect square: see the Lemma and the Theorem on p. 506 of |2Pj : 
this is a more serious obstacle. On the other hand, given an initial solution (x , y , zq) of 
(|13|) . by that Lemma one can still generate infinitely many solutions (x n , x n+ \, zq) with the 
same value of z by setting x\ = yo and iterating the integrable recurrence 

< + z 2 + J (15) 

(with N 2 Zq > 4 for an infinity of integer solutions) or, equivalently, the linear recurrence 



x n+ \ — NzQX n + x 



n-l 



0. 



(16) 



Halburd's Diophantine integrability test (THj applies to rational maps with iterates and 
parameters taking values in Q or a number field. For x G Q with x = p/q as a fraction 
in lowest terms, the height is H(x) = max{|p|, \q\}, while the logarithmic height is h(x) = 
log H (x). A rational map of this kind is defined to be Diophantine integrable if its iterates 
x n are such that h(x n ) grows no faster than a polynomial in n. To apply the test to the 
recurrence (JHJ) we shall restrict ourselves to the case when all of the iterates r n are integers, 
so that h(r n ) = log |r n |. More generally we can consider real or complex values of r n , and set 
A n = log |r n |. Now suppose that A n ~ CX n — » oo as n — > oo for real A > 1 and some C > 0. 
Then clearly r^ 1 — > and also A n+1 — A n+2 ~ CA n+1 (l — A) — >• — oo, whence T n+ i/r n+2 — > 
as n — » oo. So taking the logarithm of (JHJ) and substituting in these asymptotics gives 



A 



n+3 



2A n+2 + A n = log 



r, 



ii+i 



J 



T, 



n+2 



T, 



n+2 



0. 



(17) 



The characteristic polynomial for the linear recurrence on the left hand side of (JT7jl gives 
A 3 — 2A 2 + 1 = 0, with the largest root being the golden mean, A = (1 + y/E)/2 > 1, 
which is consistent with the original assumptions on the asymptotics of A n . Thus if integer 
sequences of this kind exist, then their logarithmic heights grow exponentially with n with 
lim n ^ 00 (log/i(r n ))/n = log((l + \/5)/2) ~ 0.4812. We proceed to show existence. 

Proposition 2 With initial data (to,ti,T2) = (1,1,1) the recurrence (OJ) gives a sequence 
of monic polynomials r n = p n {N) G Z[iV] with N = J + 3, and degp n = f n — 1 for n > 
where f n is a Fibonacci number. For real N > 2, p n+ i(N) > p n (N) > 1 holds for n > 3. 
Furthermore, for N > 2 and all n > 4 these polynomials satisfy 

(N - l)^" 1 < p n (N) < N f "-\ (18) 

That these initial data give polynomials in N follows from the earlier Proposition; equa- 
tion (HU) gives N = J + 3, and we have p 3 = N - 1, p A = N 2 - N - 1, p 5 = N 4 - 2N 3 + N-1, 
p 6 = N 7 -3N 6 +N 5 +3N 4 -2N 3 + 1 etc. From (O, d n = degp n satisfies d n+3 = d n+2 +d n+1 + l 
which implies that d n = f n — l, with f n being a Fibonacci number (/ = 0, fx = 1), and thep n 
are clearly monic. The values p n (— 1), p n (0) and p n (l) cycle with periods 4,6 and 12 respec- 
tively, while p n {2) = 1 for all n, but for iV > 2 we have J > —1 so by induction Pn +1 + J > 
and © gives p n+3 p n > p 2 n+2 so p n +z/Pn+2 > {Pn+2/Pn+i){Pn+i/Pn) > 1 as required. The 
upper bound in (JTSJl is easy to see by induction from (|T2"j) . since p n+3 < iVp n+2 p n+ i . The 
lower bound follows from another induction, by writing p n+1 = p n p n , and noting that for 
N > 2 the ratios satisfy p n > (N — l)^- 1 when n > 3, as a consequence of the inequality 
Pn+2 > Pn+iPn- Our main result is an immediate consequence of the preceding facts. 

Theorem For each integer N > 3, every triple (p n (N),p n+ i(N),p n+ 2(N)) of adjacent iter- 
ates in the sequence for n > constitutes a distinct solution triple for the Diophantine equa- 
tion (Q3J) with J = N — 3. Moreover, for all such JVgZ the logarithmic height h(p n (N)) of 
these iterates grows exponentially, and for all real N > 2 we have lim n ^ 00 (loglogp n (iV))/ri = 
log A ~ 0.4812 where A is the golden mean. 

The integrality of the terms p n (N) for iV G Z is obvious. The fact that these are distinct 
solution triples follows from the monotonicity of the sequence for N > 2, while the statements 
about asymptotic growth follow from taking logarithms of each side of (JTSJ). Similar estimates 
hold for \p n \ when iV < —2, but then p 3n (N) < 0. It is quite easy to see that log A is also 
the value of the algebraic entropy for the recurrence. 

It is hoped that the connections between Diophantine integrability, singularity confine- 
ment and the Laurent phenomenon will become clearer in the future. 
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